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Abstract 
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INTRODUCTION 



In recent years, a number of theoretical results pertaining to systems evolving far from 
thermal equilibrium have been derived; see, for instance, Refs. Jl|-§|. While not identical, 
these results bear a similar structure, and the term fluctuation theorem has come to refer 
to them collectively. Our purpose in this paper is to illustrate these and other results, 
with a highly idealized but exactly solvable model of a system far from equilibrium. 

The paper is organized as follows. In Section | we briefly review the four theoretical 
predictions which we will illustrate with our model. These are: the steady-state and transient 
fluctuation theorems, a detailed fluctuation theorem, and a nonequilibrium work relation 
for free energy differences. In Section |I| we introduce and solve our model, representing 
a particle dragged through a thermal environment by a uniformly translating harmonic 
force. In Section |TT| we use this solution to show that our model indeed satisfies the above- 
mentioned nonequilibrium results. 



Statements of the fluctuation theorem (FT) which have appeared in the literature can 
be expressed by the following equation: 



Here, a is the average rate at which entropy is generated during a given time interval 
- or segment - of duration r, for a system away from thermal equilibrium; p T (a) is the 
distribution of values of a over a statistical ensemble of such segments; and Boltzmann's 
constant ks = 1. (Note that the words "average rate" denote the time average over a 
given segment, not an average over the ensemble of segments.) Following the literature, 
we will distinguish between two versions of the FT, steady-state and transient, which differ 
in the statistical ensemble of segments considered. In the steady-state case [0, we imagine 
observing the system in a nonequilibrium steady state for an "infinite" length of time, which 
we chop up into infinitely many segments of duration r, and we compute the average entropy 
generation rate, a, for each segment; p T (a) is then the distribution of values of a over this 
ensemble of segments. In the transient case [I], by contrast, we imagine that the system of 
interest begins in an equilibrium state, but is then driven away from equilibrium, for instance 
by the sudden application of an external force. We observe the response of the system for 
a time r, starting from the moment the external perturbation is applied. Then p T (a) is the 
distribution of values of average entropy generation rate over infinitely many repetitions of 
this process. The transient FT is valid for any duration r, whereas the steady-state FT 
becomes valid as r — > oo; hence the parenthetical appearance of that limit in Eq.[TJ 

In the examples considered in the literature, the physical origin of entropy generation is 
the exchange of heat between the system and some infinite reservoir. To model the evolution 
of a system in contact with a reservoir, a variety of schemes are available, and in each case 
one must define what is meant by the "entropy generated". The FT has been derived 
or illustrated numerically for a number of such schemes, involving both deterministic and 
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stochastic equations of motion. Of particular relevance for the present paper is the work of 
Kurchan ||, who showed that the FT is valid for Langevin processes. 

Apart from the usual (steady-state and transient) statements of the FT, a detailed fluc- 
tuation theorem (DFT) was derived in Ref . f7| . This result pertains to a finite time interval 
r, and has a structure similar to that of the usual FT, but in addition exhibits a dependence 
on the initial and final microstates of the system. Specifically, consider a process Il + during 
which a system of interest evolves in contact with a heat reservoir^], while - possibly - some 
work parameter A (e.g. an applied field) is being manipulated externally. Let A + (t) be the 
externally imposed time-dependence of this parameter, from t — to t — r, and let As 
denote the entropy produced during a particular realization of this process. Now consider 
the "reverse" process, II - , defined exactly as but with the time-dependence of the work 
parameter reversed: X~(t) = A + (r — t). For the forward process IT 1 ", let P + (zf, As|zj) denote 
the joint probability that the entropy produced over the interval of observation will be As, 
and the final microstate of the system will be Zf, given an initial microstate z^; and let 
P_(zj, As|zj) denote the same for the reverse process, II - . Then the DFT states that these 
joint, conditional probability distributions satisfy the following relation: 



where the asterisk (*) denotes a reversal of momenta: (q, p)* = (q, — p). 

Finally, in Ref. M the following nonequilibrium work relation for free energy differences 
was established: 



(See also Ref. [f| for an elegant derivation of Eq.| specific to stochastic processes.) This result 
applies to a situation in which a system of interest, in contact - and initially in equilibrium 
- with a heat reservoir at temperature evolves in time as an external parameter A is 
switched at a finite rate from and initial to a final value, say, from to 1. The finite-rate 
switching of A drives the system out of equilibrium. W is the work performed during one 
realization of this process. The precise value of W will depend on the microscopic initial 
conditions of both the system of interest and the reservoir; by repeating the process infinitely 
many times (sampling initial conditions from equilibrium ensembles), we obtain a statistical 
ensemble of microscopic realizations of the process. The angular brackets then denote an 
average over this ensemble of realizations, and AF is the free energy difference between the 
equilibrium states (at temperature /3 _1 ) corresponding to the values A = and A = 1.0 EqfJ 



1 In fact, the formulation of Ref. is somewhat more general than that described here: the 
system is allowed to be in contact with numerous reservoirs, at various temperatures; furthermore, 
contact between the system and any of the assorted reservoirs can be externally established and 
broken during the course of the process. 

2 The validity of Eqj3| does not require that the system end in the equilibrium state corresponding 
to A = 1, only that such an equilibrium state exists, and that the final value of the external 
parameter is 1. 




exp(AS/k B ), 
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is valid regardless of how slowly or quickly A is switched from to 1, hence applies even if 
the system is driven far from equilibrium by a violent variation of the external parameter. 

Eqs.[I], |2], and [3] represent general theoretical predictions applicable to systems far from 
thermal equilibrium. We now introduce and solve a simple model which illustrates these 
results. 

II. THE MODEL 

In this section we introduce our (one-dimensional) model of a particle dragged through 
a thermal medium, and then we solve it exactly: given the particle's initial location, we 
compute the joint probability distribution of the final location and the net external work 
performed on the particle, after an arbitrary time of evolution. In the following section we 
use this result to confirm both the steady-state and transient fluctuation theorems, as well 
as the detailed fluctuation theorem and the nonequilibrium work relation for free energy 
differences. 

In our model, the particle under consideration obeys Langevin dynamics. This model can 
therefore by viewed (with qualifications, see Section |11 A[ ) as a special case of the situation 
considered by Kurchan |§, who showed that the FT is satisfied for this class of stochastic 
dynamics. In essence, we are illustrating Kurchan's results with a specific model. However, 
the definition of entropy generated which we choose differs from that of Ref. 0, and so 
we are obliged to express the steady-state and transient FT's differently, in terms of power 
delivered rather than entropy generation rate. We stress that this difference (between the 
formulation of the FT given below, and that in Kurchan's work) is only one of terminology, 
not substance. 

Consider the following situation. A particle, in contact with a thermal medium at tem- 
perature /3~ l , is pulled through that medium by a time-dependent external harmonic force. 
Assuming a single degree of freedom, let x denote the location of the particle, and let 

k 

U(x,t) = -(x-ut) 2 (4) 

be the moving potential well which drags the particle. We can picture the particle as being 
attached to a spring, the other end of which moves with a constant speed u. Assume 
furthermore that the thermal forces can be modeled as the sum of linear friction and white 
noise, and that the motion of the particle is overdamped. Then the equation of motion for 
the position of the particle is: 

k 

x = — (x — ut) + v, (5) 

7 

where x = dx/dt, 7 is the coefficient of friction, and v{t) represents delta-correlated white 
noise with variance 2/(3^ (as mandated by the fluctuation-dissipation relation): 

(v{t 1 )v(t 2 )) = (2/(3 1 )-5(t 2 -t 1 ). (6) 

Imagine that we observe the evolution of such a particle over a time interval from t = 
to t — t, and from the observed trajectory we compute the total work W performed by the 



4 



external potential over this interval. Then the central result of this section will be an answer 
to the following question. Given an initial location xq, what is the joint probability 
distribution for the final location and value of work performed (x, W)l 

To answer this question, we first introduce a "work accumulated" function, w(t), which 
gives the work performed on the particle up to time t; hence, W = w(r). This function 



satisfies 10 



dU 

w(t) = —(x(t),t) = -uk(x(t) -ut), (7) 

along with the initial condition w(0) = 0. 

It will furthermore prove convenient to specify the location of the particle by a variable 
y = x — ut (i.e. in the reference frame of the moving well), rather than x. Under this change 
of variables, Eqs.|| and |7] become: 

y = y — u + v (8a) 

w = —uky. (8b) 

Now imagine a statistical ensemble of such particles, represented by an evolving probability 
distribution f(y,w,t). Eq|| then translates into the Fokker-Planck equation 

= tL( y f) + U ^- + uky^- + ^-^1 (9) 
dt 7 dy dy dw flj dy 2 

What we now want is an expression for f(y, w, t\y ), by which we mean the solution to Eq.|| 
satisfying the initial conditions 

f(y,w,0\y ) = 5(y-y )5(w). (10) 

The function f(y,w,t\y ) is the joint probability distribution for achieving a location y and 
a value of work accumulated w, at time t, given y at time 0. By evaluating this solution 
at t — t, and making the change of variables from y back to x, we have the answer to the 
question posed earlier in boldface. 

To solve for f(y,w,t\y ), we first note that Eq.| has the following property: if at one 
instant in time the distribution happens to be Gaussian, then it will remain Gaussian for all 
subsequent times. This follows from the fact that the drift and diffusion coefficients in Eq.|9| 
are either constant or linear in y and w. Now, a normalized Gaussian distribution f G (y,w) 
is uniquely defined by the values of the following moments: 





1 f G (y,w) 


y 


(11a) 


w = 


f f G (y,w) 


w 


(lib) 




1 f G (y,w) 


(y 2 - f) 


(11c) 




f f G (y,w) 


(w 2 — w 2 ) 


(lid) 


c yw = j 


f f G (y,w) 


(yw - yw), 
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where the integrals are over (y, u>)-space. An explicit expression for f G (y,w) in terms of 
these moments is: 



where 



f G (y,w) = X£ exp(-z T Cz/2), (12a) 

Z7T 




Col -Cc 



-CCy W C(Ty 




(12b) 



C = (cfyCr^ — c^ w )~ l = det C, and z T denotes the transpose of z. 

The evolution of a time- dependent Gaussian distribution f G (y,w,t) is thus uniquely 
specified by the evolution of the moments y, • ■ ■ , c yw . Given a distribution evolving under 



Eqfl we get from Eq.ll the following set of coupled equations of motion for these moments: 



d A k „ d „ 

t:?/ = — V - u , — w = -uky (13a) 
at 7 at 

Jt j2 y = ~ 1 ^ o2 y~^d ' ^ c ^ = -^ c ^ + ^ a 2) > ^ = -2^ Cj/w . (i3b) 



The distribution f(y,w,t\yo) will thus be a Gaussian whose moments evolve according 
to Eq.[L3], and satisfy the initial conditions: y = y , w = <jy = c yw = <r^ = 0. The solution, 
as easily verified by substitution, is: 

y(t\y ) = -l + (y + l)e- kt ^ (14a) 
w(t\y ) = uklt + U1 {y Q + l)(e~ kt ^ - 1) (14b) 

° 2 y(t\y ) = -1(1 - e~ 2 ^) (14c) 



c yw (t\yo) = -^(e- kt/ ^l) 2 (14d) 
<£(*lw>) = - e " 2W7 + ^' ktH - 3 )' ( 14e ) 



where 



The combination of Eqs.|l4| and p72| gives the solution for f(y,w,t\y ) for all times £ > 
(and for all values of y, w, and yo). The explicit expression given in the Appendix. Note 
that, by projecting out either of the two independent variables y and w, we can obtain the 
marginal probability distributions for the other: 

p(y, t\y ) = fdw f{y, w, t\y ) = J— exp[-(y - y) 2 /2al] (16a) 



rj(w,t\y ) = f dy f(y,w,t\y ) = j= — 



■.exp[-(w-w) 2 /2a 2 w ]. (16b) 
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It is instructive to consider the limit of asymptotically long times, > 1. In this 
limit, 



y{t\y ) -> -/ (17a) 

w(t\y ) ^ uklt + 0(1) (17b) 

a 2 y (t\y ) - 1/pk (17c) 

c^(t|yo) -> -wy/pk (17d) 

^(t|y )-^2 U V//9+O(l), (17e) 



where the "order unity" corrections to w(t|?/o) and <r^(t|?/o) represent terms which converge 
to a constant as t — > 00. We see that, for any y , the distribution of positions settles 
into a steady-state Gaussian of variance l/[3k, centered at a displacement — / from the 
minimum of the confining potential U (Eqs. |17a| , |1 7c|) . Note that a canonical (equilibrium) 
distribution of positions would have the same variance, but centered at the minimum. Hence, 
I represents the extent to which the steady-state distribution of positions "lags behind" the 
instantaneous equilibrium distribution. This lag could have been predicted with an educated 
guess, by ignoring fluctuations and simply balancing the frictional and harmonic forces, — 'ju 
and — ky. One can similarly understand the leading behavior of w(t): if the position of the 
particle is (on average) a distance I behind the instantaneous minimum of the well, then the 
harmonic spring is pulling the particle with a force +kl, at a velocity u, hence delivering a 
power ukl. 



A. Energy conservation and entropy generation 

In our model, energy conservation translates into the following balance equation between 
the external work W performed on the particle, the heat Q absorbed from the thermal 
surroundings, and the net change in the internal energy of the particle [10]: 



W + Q = AU, (18) 
where AU = U(x(t),t) - u(x(0),0), W = w(f), and 

Q = £ x(t)^-(x(t),t)dt (19) 

In addition to these quantities (W, Q, and AU), we need a microscopic definition of the 
entropy generated over one realization of the process. There is necessarily some arbitrariness 
in such a definition, but we will use the following one: 

AS = -PQ. (20) 

Thus, we identify the entropy generated with the amount of heat dumped into the envi- 
ronment (-Q), divided by the temperature. This definition - motivated by macroscopic 
thermodynamics (see e.g. Ref. [0]) - differs from that of Kurchan pj, who identifies the 
entropy production rate with the external power delivered to the system, divided by the 
reservoir temperature. Both definitions seem to be reasonable, but we will use Eq.|20|. 
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Quite apart from the definition of entropy generation, our set-up differs from Kurchan's 
in the way in which external work is defined; see Section 2.2 of Ref. ||. Nevertheless, 
the general analytical approach taken by Kurchan applies to our set-up as well. Hence, 
any results derived in Ref. regarding average power delivery (external work performed, 
divided by the duration of the time interval), extend to our situation. 

III. ILLUSTRATION OF FAR-FROM-EQUILIBRIUM PREDICTIONS 

We now use the results of the previous section to show that Eqs.[l|-|3| are obeyed by our 
model. 

A. Transient and steady-state fluctuation theorems 

We begin with the transient and the steady-state fluctuation theorems. Ordinarily the 
FT is expressed in terms of the average entropy production rate (Eq.[l]). By contrast, the 
relations which we will show to be satisfied by our model (Eq.|2"2"|) are expressed in terms of 
the average power delivered to our particle as it is dragged through its thermal environment. 
This difference - as mentioned earlier - is a consequence of our choice of definition of 
average entropy generation rate, a = AS/r = —f3Q/r. (Kurchan, by contrast, defines 
entropy generation in terms of power delivered, (XKurchan = PW/r. ||) In the calculations to 
follow, the reader should bear in mind that the transient and steady-state relations which we 
establish are essentially those of Kurchan; only our definition of entropy production prevents 
us from presenting them as such. 

(It might be interesting to investigate whether or not EqlJ remains valid under our 
definition of entropy production, a = —f3Q/r. This would be easy to check with numerical 
simulations, but in the "exactly solvable" spirit of restricting ourselves to analytical results, 
we have not pursued this question.) 

For a time interval of duration r, let 

X = W/r (21) 

denote the time-averaged rate at which work is performed on the particle - i.e. the average 
power delivered - by the moving harmonic potential. Now imagine that we observe the 
particle for an "infinitely" long time as it evolves in the steady state; we divide this time 
of observation into infinitely many segments of duration r; we compute the average power 
delivered, X, over each segment; and we construct the statistical distribution of these values, 
p^(X). We will obtain an explicit expression for pf(X) for our model, and will show that it 
satisfies the following steady-state FT: 

lim -ln ^±^ =/3X. (22a) 

Now imagine instead that we begin with the particle in thermal equilibrium, and then we 
drag it for a time r with the harmonic confining potential. (That is, for t < the potential 
is stationary, with the particle in equilibrium; then, between t — and t — r the potential 



S 



is moved rightward with velocity u.) Again denning X to be the average power delivered 
over the interval < t < r, let p^(X) denote the statistical distribution of values of X, over 
infinitely many repetitions of this process, always starting from equilibrium. [] We will solve 
for p^'(X) and show that it obeys the following transient FT: 



1, P?( 

r p^{ 



+X) 



-x) 



px, 



(22b) 



whose validity does not require the limit r — > oo. 

Let us begin by considering the steady-state case. Since the right sides of Eqs. |17a| and 



17c| are independent of yo, an arbitrary initial distribution of particle positions p(y, 0) will 
settle into a Gaussian: 



lim p(y, t) 



t— too 




-f3k{y + lf/2]=p s {y), 



(23) 



which defines the instantaneous distribution of particle positions in the steady state. Then 
pf{X) - the distribution of values of average power delivered, X = W/r, over time intervals 
of duration r sampled during the steady state - can be constructed by folding together 



P (l/o) and T](w, r\y ) (Eq.fB§): 



pf{X) = J dw5(X -w/t) J dy p s (y )ri(w,r\y ). 



(24) 



Here, the integral over dy produces the distribution of values of work, after time r, given 
initial conditions sampled from the steady state. The integral over dw converts that dis- 
tribution of values of w into one of values of X. Examining Eqs. |14b| , |14e , |16b| , and we 
see that, in the product p s (y )r](w, r\y ), the variable y appears only in powers up to the 
quadratic inside an exponent; hence this product is a Gaussian in y , and the integral can 
be carried out explicitly. Without going through the (modestly tedious) details, we present 
the result: 



P S T (X) 



1 



2na 2 x 



cxp 



{x - xy 

2a\ 



(25) 



where 



X = ukl 



a 



x 



2fiX 



1 + ^{e~ kTh 
kr 



1)- 



(26) 



(X represents the average instantaneous power delivered to the particle, in the steady state; 
see the comments following Eq.[l7|.) From this we obtain an explicit expression for the left 
side of Eq. |22a| : 



The superscript C indicates that the particle's initial conditions are sampled from a canonical 
ensemble. 
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1 It, 



1 2XX (3X 



a 



x 



(27) 



Since lim^oo p = 1, we conclude that the FT for power delivered in the steady state (Eq. |22a| ) 
is indeed satisfied for our model. Note, however, that the limit r — > oo is necessary. 

In the case of the transient FT, pertaining to a system driven away from an initial state 
of canonical equilibrium, we can construct p^'(X) in the same way as pf(X), only now we 
fold rj in with a canonical distribution, p (y) oc exp(— (3ky 2 /2), rather than the steady-state 
distribution: 



p?(X)= I dw6{X-w/r) I dy p c \y )ri(w,T\y ) 



2na 2 x 



exp 



(X - pXf 
2a\ 



(28) 
(29) 



where a x , X, and p(r) are exactly as above. The only difference between the steady-state 
and the transient distribution of values of X is, we see, the factor p appearing inside the 
exponent in the latter. This small difference has the effect that the FT for power delivered 
in the transient case is satisfied for all positive values of r: 



1 2pXX 



T pC(-X) 



a 



px. 



(30) 



x 



B. Detailed fluctuation theorem 

We now show that our model satisfies the detailed fluctuation theorem (DFT). As men- 
tioned in Section [I], the DFT is stated in terms of two processes, Il + ( "forward" ) and IT - 
("reverse"), related by time-reversal. In the present context, we take Il + to be the pro- 
cess studied above: a particle is dragged through a thermal medium by a time- dependent 
potential well, 

U+(x,t) = -{x-utf. (31) 

The reverse process, II - , is then obtained by moving the well in the opposite direction: 

k 

U~(x,t) = -(x + ut-ur) 2 . (32) 

Formally, we can think of the minimum of the potential as being given by XAx, where 
Ax = ur is a constant, and A is an externally controlled parameter. During the process II + , 
A is changed uniformly from to 1; during II - , from 1 to 0. 

Again taking y to be the displacement of the particle relative to the instantaneous min- 
imum of the potential, we define f + (y,w,t\y ) and f_(y,w,t\y ) to be the joint probability 
distributions for attaining (y,w) at time t, given y at time 0, for the two processes. We 
then solve for these two distributions exactly as we solved for f(y, w, t\yo) in Section [II]. The 
solutions are: 
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f+(v, w, t\y ) = f(y, w, t\y ) (33) 
f-(y,w,t\y ) = f(y,w,t\y ) u ->- u , (34) 

where / is just the solution obtained in Sec.||. Thus, the solution for the forward process 
is identical to the solution of Sec.fT] (as it must be, since Il + is exactly the process studied 
there!), whereas the solution for the reverse process is obtained from / by replacing u by 
— it everywhere, including in the definition of I (hence, I — > —I). This replacement is easily 
understood: in terms of the variable y, the process II - is no different than that obtained by 
starting with the minimum of the potential at x = 0, and moving it with velocity — u for a 
time t. 

Let us now put these solutions to good use. 

The DFT, in the context of this problem, claims the following: 

P+(yB,+AS\y A ) 

— — j — r = exp AS, 35 

P-{y A ,-AS\y B ) 

where P±(yf,As\yi) denote the joint probability distributions of finding the particle at a 
final point yj, and a value of entropy generated As, given an initial location y^, for the 
forward and reverse processes.^ Now consider a single realization of either process. Energy 



conservation, combined with our definition of entropy produced (Eqs.[0| and |20), give us the 
following relation between y^ yf, As, and the work w performed: 

w - /T'As = k -(y) - y\) = AU( yi ,y f ). (36) 
We can then re-express P± in terms of f±: 

P ± (y f ,As\ yi ) = J dwf ± (y f ,w,T\ yi )5(As + l3AU-l3w) (37) 

= r 1 f±(y f , \{y) - y- ) + r 1 ^, r\ yi ), (38) 

from which we obtain the following explicit expressions for the numerator and denominator 
in Eq.|BJ: 

P+(y B , +AS\y A ) = p-'Uiys, r\y A ) (39) 
P4y A ,-AS\y B )=/3- 1 f-{y Ai -n i T\y B ) } (40) 



where 



Sl = ^{y B -y A )+(3- l AS. (41) 



Eq.|35| is thus equivalent to the following relation 

exp 



f-(y A , -Q,r\y B ) 



2 



where, if Eq.^ is to be valid for all real AS, then Eq.^2| must hold for all real Q. Using the 
expression for / given in the Appendix, one can verify that, indeed, Eq.|42| is valid for all 
values of Q, hence the DFT is satisfied by our model. 



4 Strictly speaking, these should be defined in terms of absolute locations x A and x B , but the 
change of variables to relative displacements y A and y B is immediate. 
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C. Generalization 



Before proceeding to the nonequilibrium work relation, we point out that our model is 
easily generalized to include an additional uniform, constant external force. Namely, suppose 
we modify the time-dependent potentials U ± (x,t) for the forward and reverse processes, as 
follows: 

U ± {x,t)-^U ± {x,t) + ax , a > 0. (43) 

This corresponds to subjecting the particle to an additional leftward-pushing force, of mag- 
nitude a. Thus, assuming u > 0, the particle is dragged "up" the potential energy slope ax 
during and "down" the slope during n - . 

The solution in this case, for the forward process n + , is the same as that in the previous 
section, except that I is everywhere replaced by 

l> = J+ « = o + 2! ,_ (44) 

Thus, in the steady state, the average position of the particle is displaced by an amount 
a/k to the left, relative to the case with no additional force (Eq. |17a| ). This implies that 
additional work is performed on the particle at an average rate ua (Eq. |17b| ); this is simply 
the average rate at which we drag the particle up the slope. 

For the reverse process, the solution is obtained (as in the previous Section) by the 
further replacement u — ► —u, including in the definition of l a . 

It is straightforward to show that, with these replacements, the DFT remains valid. 



D. Nonequilibrium work relation for free energy differences 

Let us finally use the results derived above to show explicitly that Eq.|| is satisfied by our 
model. We will consider the process defined by U + (x, t) in the previous section (Eq.|4"3"|). As 
before, let Ax = ut be a fixed distance, and let AAx define the minimum of the confining 
potential, so that we move that minimum from to Ax by changing A from to 1: 

k 

U x (x) = -{x - XAx) 2 + ax (45) 
U + {x,t) = U m {x), (46) 

where X(t) = tjr. For any fixed value of A and fixed temperature there exists an 
equilibrium state of the system, defined microscopically by a canonical distribution. The 
associated free energy F\ is then defined in terms of the logarithm of the corresponding 
partition function: 

F x = -/T 1 In J dx e -^ (x) (47) 
a 2 1 , Pk . . 

= aXAx -Tk + 2fJ hL ^ (48) 

Hence the free energy difference is simply 
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AF = F x - F = aAx. 



(49) 



Physically, this is the work required to reversibly (i.e. infinitely slowly) change A from to 
1, at constant temperature. 

Now consider a statistical ensemble of realizations of our finite-time, irreversible process, 
X(t) = tjr, with initial conditions sampled from the canonical ensemble corresponding to 
A = 0. The associated distribution of values of work, W — w(r), can be written as: 

V (W) = J dy pC(y ) V (W,T\y ), (50) 

where p^iuo) is the canonical distribution of initial conditions (for a^O), and r](W, r\y ) is 
the distribution of work values at time r, given initial condition y . We can use the results of 
Section [Tl], with the substitutions / — > l a and u = Ax/t, to compute this integral explicitly. 
Once again skipping the algebra, we present the result: rj(W) is a Gaussian distribution of 
mean and variance 

(W) = aAx + ^(Ax) 2 /t (51) 
a 2 w = 2w(Ax) 2 /Pt, (52) 

with /i(r) as defined by Eq.^6[ (Note that in the reversible limit - i.e. r — > oo with Ax fixed 

- we get (W) = aAx and a\y = 0, hence W = AF for every realization, in agreement with 
the remarks following Eq.^.) Thus, for any positive value of r, r](W) is a Gaussian whose 
mean and variance are related by 

(W) = AF + (3a 2 w /2. (53) 

This implies - as can be verified by direct evaluation of / dWt](W) exp(— (3W) - that the 
nonequilibrium work relation for free energy differences, Eq.[3] above, is satisfied. 

Eq.[53] has the structure of the usual fluctuation- dissipation theorem (FDT) for linear 
response, but differs from the latter in the fact that the fluctuations - represented by a^y 

- are not obtained from the equilibrium fluctuations of the particle, but rather truly reflect 
behavior far from equilibrium. Indeed, the presence of //(t) in Eq.[5| is evidence of non-linear 
response in our model: for fixed Ax, the average dissipated work, (W) — AF, is not simply 
inversely proportional to r, as would be the case for linear response. 
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APPENDIX 



Here we give an explicit expression for the function f(y,w, t\y ) introduced in Section O. 
This solution is essentially the combination of Eqs.O and 111: 



f(y,w,t\y ) = ^exp(AC/2), (54) 

Z7T 



where 



B 2 k 2 

C = detC = - - £ (55) 

27u 2 z/_[27z/_ + ktu + ) v ; 

A = (l//3k){jy + v-{w - kltu) 2 + 7 2 uV(y - y )[Alv- + (3v - l)y Q + [y - 3)y] (56) 
+2^u{l 2 kutv 2 _ — \v_\2wv_ + ktuv + (y — y)] — wu 2 _(y + y) — ktu(y — z^/o) 2 )}; (57) 

and 

v = exp(— kt/j) , u± — v ± 1. (58) 
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